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Monte Carlo simulations of ordering in ferromagnetic-antiferromagnetic bilayers 
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Monte Carlo simulations have been used to study phase transitions on coupled anisotropic 
ferro/antiferromagnetic (FM/AFM) films of classical Heisenberg spins. We consider films of different 
thicknesses, with fully compensated exchange across the FM/AFM interface. We find indications of 
a phase transition on each film, occuring at different temperatures. It appears that both transition 
temperatures depend on the film thickness. 



INTRODUCTION 

Most of the attention given to ferromag- 
netic/antiferromagnetic (FM/AFM) bilayers is normally 
directed to the changes of the magnetic properties of the 
FM when it is in contact with the AFM. The commonly 
observed effects are a unidirectional shift (exchange 
bias) and a significant increase of the coercivity [lj . The 
blocking temperature, Tb, below which these effects are 
observed, is comparable to the bulk Neel temperature 
when the AFM film is thick, but can be considerably 
lower when the AFM film is thin 0, H S H @. It is 
generally believed that this reduction of the blocking 
temperature is due to finite size effects in the AFM, i.e. 
the ordering temperature of the AFM is decreased due 
to its finite thickness . In a recent neutron experiment 
on CoO/Fe304, however, van der Zaag et. al. |j| have 
found that the AFM shows signs of ordering above Tb- 
This indicates that the proximity of the FM influences 
the phase transition in the AFM in a way that cannot 
be predicted from studying free AFM films. 

In this paper we use Monte Carlo simulations to an- 
alyze the effect of the FM on the transition between 
disordered and ordered states in the AFM for various 
film thicknesses and compensated exchange across the 
FM/AFM interface. 



MODEL AND METHODS 

The system studied here consists of a multi-layered fer- 
romagnetic (FM) film coupled to an underlying multi- 
layered antiferromagnetic (AFM) film with no lattice 
mismatch at the FM/AFM interface. The Hamiltonian 
of the model is given by 

n = - j f J2 Sr ' Sr ' - ^E (^) 2 

(r,r')GFM reFM 

-J A J2 Sr-Sr'-K A J2 ( S r) 2 
(r,r')eAFM rGAFM 



-Jaf ^ S r 

(r,r')eFM/AFM 



(1) 



where S r = (S* , S% , S') is a three-dimensional classical 
Heisenberg spin of unit length, (r,r'} denotes nearest- 
neighbor pairs of spins coupled with exchange interac- 
tions J F > on the FM film, J A < on the AFM film, 
and Jaf at the FM/AFM interface, which is on a (001) 
plane. Spins on the AFM film have a uniaxial single- 
site anisotropy K A > 0, whose easy axis is along the y 
axis. In contrast, the single-site anisotropy for spins on 
the FM film has a hard-axis (Kp < 0) along the z direc- 
tion, which is perpendicular to the FM/AFM interfacial 
plane. No external magnetic field is applied on the films. 
The structure of the films is a body-centered cubic lat- 
tice, with linear sizes L x , L y , and + Lf , measured 
in terms of two-spin unit cells. and denote the 
number of unit-cell layers on the AFM and FM films, re- 
spectively. Hence, the total number of spins in the lattice 
is N = 2L x L y (Lf + Lf). We use periodic boundary con- 
ditions along the x and y directions and free boundary 
conditions along the z direction. 

The effect of the film thickness on the FM and AFM 
transition temperatures is studied for fixed interaction 
parameters. We consider Jp = J > 0, J a = — J, 
Ka = J, Kp = — 0.5J, and a compensated interface 
with Jaf = r J, where r is a random number uniformly 
sampled in the interval [—1,1]. Values of the film thick- 
ness used are = Lf = 3, 6, and 12, with several 
cross sections (L x — L y — 12, 20, 40, and 60) to analyze 
finite-size effects. 

The order parameter for the ferromagnetic transition 
is the uniform magnetization per spin m — | ^ r S r \/N, 
whereas to characterize the antiferromagnetic transition 
it is necessary to divide the BCC lattice into two sim- 
ple cubic sublattices, denoted I and II, and consider 
the staggered magnetization per spin defined as m s — 
l£ re /S r -£ re //S r |/iV. 

Our simulations were carried out using importance 
sampling Monte Carlo methods @ , with Metropolis algo- 
rithm, at fixed temperature T. Typically 3 x 10 5 Monte 
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Carlo Steps/site (MCS) were used for computing aver- 
ages after about 1 x 10 5 MCS were discarded for thermal- 
ization. Whenever not shown, error bars in the figures 
are smaller than the symbol sizes. 



RESULTS 

Figs. Jit) and (QJ)) show the uniform and staggered 
magnetizations per spin for films of different thicknesses 
and fixed cross section. At low temperature, both m and 
m s tend to 0.5, indicating ordered FM and AFM films 
(note that both quantities are computed for the entire 
lattice and thus go to 0.5 instead of 1.0 when the films 
are ordered). As the temperature increases, both m and 
m s decay to zero, indicating disordered spin configura- 
tions on both films at higher temperatures. In simu- 




FIG. 1: (a) Uniform and (b) staggered magnetizations as 
a function of temperature for several film thicknesses, with 
L x =L y = 20 and = L F Z = L z . 

lations of finite lattices, occurrence of phase transitions 
cannot be ascertained with the use of only one lattice size. 



Therefore, we have considered different cross sections for 
each film thickness. Illutrations of finite-size effects on 
the uniform and the staggered magnetizations near the 
phase transitions are presented in Figs.J30 and re- 
spectively. These figures show that as the cross section 
increases, the decay of the magnetizations m and m s to 
zero becomes sharper. Such dependence of order param- 
eters on finite lattice sizes is characteristic of real phase 
transitions. Similar finite-size effects were observed for 
the uniform and the staggered magnetizations of films 
of other thicknesses. The behavior of the uniform and 
staggered magnetizations shown in Figs.QJi) and QJ)), 
with finite-size effects described above, suggests that for 
a given film thickness there are two distinct phase tran- 
sitions, one occuring on the AFM film and one on the 
FM film. The AFM transition temperature seems to be 
slightly higher than the FM one, presumably due to the 
higher anisotropy on the AFM film. In addition, it ap- 
pears that the transition temperatures on both the FM 
and the AFM films increase with film thickness. We use 



0.20 



§ 0.15 

N 

1)0.10 

c3 



| 0.05 

'§ 

3 



0.00 

1 



~f — 
o+ 



□ o 



— 1 1 <~ 

+ L =L =12 

x y 

o L =L =20 

" x y 

a L =L =40 

x y 

□ L=L=60 

x y 



o 

□a O 



(a) 



oo 



o OO 



D □ □ 6 6 



° o O O 



2.0 



2.2 
k B T/J 



2.4 



2.6 



a 



C3 



0.20 



0.15 



I, 0.10 



T3 



bo 
M 



0.00 



1 % 1 


i 


+ 

a 

\ 

- 1 


+ L =L =12 

x y 

o L =L =20 

" x y 

a L =L =40 

x y 

□ L=L=60 

x y 


A ° 

- (b) ^ o 


n + - + . 

° o , 


i 


T ¥ ' 



2.0 



2.2 



2.4 
k B T/J 



2.6 



2.8 



FIG. 2: (a) Uniform and (b) staggered magnetizations as a 
function of temperature, for = Lf = 6 and several cross 
sections. 
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crossings of reduced fourth-order cumulants of order pa- 
rameters 0, defined as U± r = 1 — (m i )/(3(m 2 ) 2 ) and 
f7| r = 1 — (m^)/(3(m 2 ) 2 ), to locate the phase transi- 
tion temperatures on the FM and AFM films, respec- 
tively. Preliminary analyses of cumulant crossings in- 
dicate that the AFM transition for Lf = = 6 oc- 
curs at T c = (2.235 ± 0.005) J/ks, and it is consistent 
with the two-dimensional Ising universality class. For 
the same film thickness, the FM transition appears to be 
at T c — (1.95 ± 0.01)J/fcs; however, the nature of this 
transition has not been determined yet. 

Figs.(I3K) an d (03) show the x, y, and z components 
of the uniform magnetization per site as a function of 
film layer, for Lf — Lf = 6, at T = 0.6J//cs and 
T = 2.0J/k B , respectively. While Fig.lJSt) illustrates 
the behavior of the layer magnetizations when both films 
are in the respective ordered states, Fig. (O^) corresponds 
to a temperature above the FM transition and below the 
AFM one. In our notation, layers 1 to 12 belong to the 
AFM film, whereas layers 13 to 24 comprise the FM film. 
Layer magnetizations on the AFM film are almost en- 
tirely along the y-axis, which is the easy axis for the 
uniaxial anisotropy on this film. Because of the BCC 
lattice structure, spins on consecutive layers belong to 
different sublattices; hence consecutive layer magnetiza- 
tions on the AFM film have different sign. In contrast, 
layer magnetizations on the FM film in the ordered state 
are essentially parallel, with a large component along the 
x-axis. In the absence of the AFM film, spins on the 
FM film have global rotation symmetry on the x-y plane, 
which is the easy plane for spins on this film. However, 
coupling to the AFM film causes spins on the FM lay- 
ers to orient preferentially in a direction perpendicular 
to the AFM easy axis. The slightly lower values of the 
magnetizations for layers 1 and 24 in Fig.J2tO are due to 
the free boundary conditions used along the z direction. 
Similar boundary effects are seen in Fig. (J2t>) ■ 
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FIG. 3: Components of the uniform magnetization as a func- 
tion of lattice layer, for — if = 6, L x = L v = 20 and (a) 
T = 0.6J/fcs and (b) T = 2.0J/fc B . The solid lines are guides 
to the eyes. 



CONCLUSIONS 

We have used extensive Monte Carlo simulations to 
study a system of coupled FM/AFM films, with com- 
pensated exchange across the interface. For the values of 
anisotropies and interaction parameters considered, it ap- 
pears that the FM and the AFM phase transitions occur 
at different temperatures. The latter transition is con- 
sistent with the two-dimensional Ising universality class. 
Our preliminary simulations suggest that both phase 
transition temperatures increase with the film thickness. 
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